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I.  INTRODUCTION 


-  Appropriate 
Mean  Radius  at  x 


I 

Particle  Surface 


Figure  1 .  The  Pore  Tree 


To  properly  describe  coupled  chemical 
reactions  and  gaseous  diffusion  in  porous 
sorbent  and  catalyst  grains,  the  "pore  tree"  was 
introduced  by  Simons  and  Finson  (1979)  and 
Simons  (1982).  The  pore  tree  represents  an 
isolated  sub-structure,  allowing  diffusion  into 
and  out  of  the  porous  media  without  permitting 
transport  through  the  media.  This  pore  structure 
w^  developed  via  analogy  to  the  kinetic  theory 
of  gases  with  the  pore  length  analogous  to  the 
mean  free  path.  Under  the  assumption  that  the 
pore  aspect  ratio  ^re  length  to  radius)  is  a 
constant,  a  pore  size  distribution  was  obtained 
that  has  been  confirmed  for  coal,  coal  char, 
sorbents,  catalysts  and  kidney  stones  from  both 
men  and  women.  The  pore  tree  was  statistically 
derived  from  the  pore  size  distribution  and 
allows  the  orderly  migration  of  a  reactant  gas 
from  the  large  pores  to  the  small  pores  (Fig  1) 

u  spatially  dependent  transport/reaction  equations  are  solved  for  a 
^gle  ^re  tree  and  then  the  total  contribution  of  all  trees  (of  all  sizes)  in  thTsy^^is 
btamed  by  summmg  the  contribution  of  each  tree  that  reaches  the  exterior  of  Ae  system 

.'s  h“”’  transpon  ^JLdon  for  a”  m»le 

”bulk”  f  ^  poi'es  at  a  fixed  point  in  space  before  integrating  spatially  The 

phcit  functions  of  pore  size.  One  example  of  this  implicit  pore  size  dependence  is  that  nf 
e  eterogeneous  reactions  within  porous  catalysts  and  sorbents  for  which  the  pore  tree 

s«„a„3po„  model  was  developed.  A  seld  example  is  hm  of  «  d^fSsfon  and 
remediation  reactions  in  the  immobile  region  of  soil.  ainusion  and 

In  order  to  describe  the  subsurface  transport  of  gas  and  water  in  soil  the  disn^rcmn  «f 
con  mmMts,  and  m-situ  remediation  of  contaminated  sites,  the  pore  tree  is’extended  herein 
o  simulate  ^rme^U,ty  and  bulk  transport.  The  random  nature  ^the  pore  s^cmre  wWch 
fo^the  basis  of  d«  sutisdcal  derivation  of  the  pore  tree,  is  applied  fo  p<^"ii  “  d 

the  "branche«^^thT^^'^^  structure  is  obtained  via  a  statistical  determination  of 

l^ge  scalffmobi^iT^‘'T"'°"  to  several  trees  to  allow  convection  and  diffrision  through  the 
ge  ^e  (mobile)  structure  m  addition  to  diffusion  and  coupled  reactions  in  the  small 

PP**  "’“‘P' 'O  ex;L  tJ!;  m™e„, 

errors  in  the  permeabihty  of  soil  due  to  the  measurement  scale  size  fShouse  et  ai  i  ooa^  a 

“  S — s^u  ^aT  r 

et.  al.,  1994)  as  a  ftmction  of  a  saturation  scale  size.  The  theory  also  depicts  a  permeable 

s'3^1  d'ffi  ^  Ih®  flow  beween  the  Llk  convection  and 

e  small  scale  diffusion.  The  model  provides  an  analytic  description  of  a  oore  sfnirtum  w 

so  upon  which  transport  and  coupled  chemical  reactions  may  be  accurately  superimposed 
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II.  ISOLATED  PORE  TREE:  THE  STRUCTURE 

Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  1^  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  W  and  is  given  by 


r„=2a0''’/3Ar„  (1) 

where  0  is  the  total  porosity  of  the  particle  and  is  a  constant  of  integration, 
approximately  equal  to  five,  which  relates  the  pore  length  to  its  radius 


The  radius  of  the  smallest  pore  is  denoted  by  r„3^  and  is  given  by 


(2) 


r^=20/pp^5^  (3) 

where  is  the  density  of  the  solid  matrix,  Sp  is  the  specific  internal  surface  area  (several 
hundred  mVg),  and 


P  =bi(r^/r^) 


(4) 


The  particle  contains  a  continuous  distribution  of  pore  sizes  from  r„,i„  to  .  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  -l-drp  is  denoted  by  g{r^)Adr^.  The  pore  distribution  function  g(r^)is  given  by 


S(r,)=9/2.tpr;  (5) 

where  g(r^)  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  27trp^.  Hence,  the  porosity  is  the 
27trp^  moment  of  5(r^)  and  the  internal  surface  area  is  the  47trp  moment  of  g(r^).  The 
expression  for  ^(r^)  was  derived  (Simons  and  Finson,  1979)  from  statistical  arguments  and 
has  been  validated  through  extensive  comparison  of  the  predicted  volume  and  surface  area 
distributions  with  mercury  intrusion  data  (Stacy  and  Walker,  1972).  This  has  been 
accomplished  for  coal,  char  derived  from  that  coal  (Kothandaraman  et.al.,  1984),  sorbents, 
catalysts  and  even  kidney  stones  from  both  men  (low  porosity  oxalate)  and  women  (high 
porosity  phosphate). 
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A  charactenstic  feature  of  the  distribution  depicts  that  the  pore  volume  between 

r  and  increases  linearly  with  the  natural  log  of  r..  It  is  the  functional  form  of  this 
relationship, 


Pore  Volume  «  j r^g(rp)dr^  «inr^ 

^miw 


(6) 


that  depicts  the  inverse  cubic  dependence  of  gir^)  on  r^.  A  linear  display  of  mercury 
intrusion  volume  vs.  Infr^)  always  infers  a  distribution. 

The  number  of  pores  within  the  bulk  volume  V  whose  pore  radius  is  between  rp  and 
Tp+drp  may  be  defined  by  Vf(r^)dr^.  The  pore  volume  is  expressed  as  the  jrrp%  moment  of 
f(r^)  and  the  internal  surface  area  is  the  27trplp  moment  of  The  pore  size  distribution 

functions  (  /(r^)  and  g(r^)  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  f(r^)  is  related  to  g(r^)  by 

8(rj,)=f(r^)l^l2  (7) 

^nation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 

The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
IS  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 
functions.  The  analysis  suggests  that  Ip ,  g(r^)  and  f(r^)  are  proportional  to  rp  ,  1 V  and 
l/Tp  respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e. 
the  total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  The  expression  for 
/(r^)  is  given  by 


f(fp)  = 


(8) 


where  the  constants  were  defined  above. 

,  .  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 

that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspect  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspect  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  r„,^  to  imply  (via  KJ  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g(r^)4na^dr^  where  f(r,)  is  functionally  identical  to  g(r^).  Each  trunk  of  radius  r,  is 
associated  with  a  specific  tree-like  structure.  Let  N,  be  defined  as  the  branch  distribution 
function  where  N,drp  is  the  number  of  pores  of  radius  rp  (within  size  range  drp)  in  a  tree 
whose  trunk  radius  is  r,  ,  The  total  number  of  pores  of  radius  rp  in  a  sphere  of  radius  a 

may  be  expressed  as  4l3na^f{r^)dr^  or,  as  the  sum  of  all  pores  of  radius  rp  contained 

within  every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the 
trunk  of  a  tree.  Hence, 


=  f  Njgir^)4na^dr^  +  4na^g(r^) 


(9) 


where  g(r^)  is  the  number  of  tree  trunks  per  unit  external  area  of  the  porous  sample  and  only 
those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r„^,  ^(r^)and  /(r^),  Eq.(9)  is  identically  satisfied  by 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  S,(rJ  and  VtCrJ, 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


(11) 


Ip 


(12) 


Using  Eq.(lO)  for  Nt,  St(rJ  and  V,(rJ  become 
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(13) 


S,(r,)  =2Tir,/^ 


V^miny 


(1-6) 


Vt(r^)  1+ln 


(14) 


where  the  (1-6)  term  in  has  been  included  to  account  for  pore  combination  (Simons,  1979). 

The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  St  and  V,  reduce  to  those  appropriate  to  a  single  cylindrical  pore  in 
the  Umit  of  rt  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  S,(r,)  and  V,(rJ  over 
®11  5  (^,)  recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 


Each  trunk  of  radius  rj  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at 
location  x  in  a  tree  of  trunk  radius  r,.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


n(x)=rflrj(x)  (15) 

and  the  coordinate  x  is  related  to  rp  by 

dr^ldx= -r^n,  (16) 


The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons,  1983b  &  1984) 
and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al.,  1986).  It  was  also 
used  to  successfully  describe  sulfur  sorption  (SO2  and  HjS)  by  porous  calcine  (CaO)  in  the 
limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984)  and  was  later 
extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Carman,  1986;  Simons  et  al.,  1987; 
Simons,  1988;  Simons  et  al.,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 
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in.  INTERCONNECTIVITY 

The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 
distribution  of  the  permeability  is  to  determine  the  distribution  function  G,(r,,rprfr  which 
represents  the  number  of  pores  of  radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  r,.  Consider  an  infinite 
homogeneous  isotropic  porous  media  and  isolate  a  spherical  volume  of  that  media  denoted  by 
the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  2.  The  total  number  of  pores  of  radius  r 
(within  size  range  drp)  intersecting  plane  AA  of  area  Tta^  has  previously  been  defined  by 

•  The  pores  in  plane  AA  in  this  size  range  may  also  be  determined  by 

integrating  Gf(r^,r^na^dr^  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the 
porous  sample.  Hence  it  follows  that 


8(rp)Tta^dr^  =  J [G^(r,,r^%a^dr^]  gir^)4%a^dr^  (17) 

where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 

A  solution  to  Eq.  (17)  for  G^(r^,r^  will  not  necessarily  be  unique.  Physical  arguments 

will  help  determine  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since 

Nt  represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore 

intersecting  a  plane  is  proportional  to  its  length,  it  follows  that  G,(r^,rp  should  be 


Figure  2.  Spherical  Volume  of  a  Porous  Media 
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proportional  to  the  product  of  N.  and  ,  i.e.,  proportional  to  r.^/r/.  Eq.  (17)  is  identically 
satisfied  by  a  function  which  differs  from  x^lr^  by  In(rp). 


(18) 


4na^r/ln(r, 


Ir  ) 

mai' 


Note  that  ln(rp/r^)  introduces  an  integrable  singularity  at  rp=r^  such  that  is 

finite  at  rp=r^  .  Hence,  there  is  one  and  only  one  largest  pore  for  each  reference  sphere. 

The  probability  of  trees  sharing  common  branches,  i.e.,  the  inter^nnectivity  of  the 

pore  structure  is  described  in  Fig.  3.  We  seek  the  distribution  function  /(r^)rfr^  which 
represents  the  number  of  pores  (within  size  range  dr,  about  rp)  per  unit  area  of  plane  AA  that 
are  connected  to  both  sides  of  the  pore  structure  through  pores  at  least  as  large  as  r  is 
defined  as  the  area  within  plane  AA  that  is  open  to  one  side  of  the  porous  media  through  all 
trees  of  size  r/  (through  all  pores  of  size  rp'  that  are  at  least  as  large  as  rp).  Subsequen  y, 

A  G(r  rjdr  represents  the  number  of  pores  of  size  rp  (within  size  range  drp)  per  unit  area 

of  plln^AA  that  are  contained  in  a  tree  of  size  range  dr,  about  r,  and  are  al^ 

the  opposite  side  of  the  porous  media  through  all  trees  denoted  by  r, .  It  follows  that  the 

distribution  function  for  interconnected  pores  in  plane  AA  may  be  obtamed  by  integrating 

enough  to  contain  a  pore  of  size  Tp  .  Hence, 


Jir^)%a^dr^  =  J  [AgG,(r,,rp)dr^  f(r,)27ra^dr, 


(19) 


Figure  3.  Interconnectivity  of  a  Porous  Media 
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From  the  above  definition  of  ,  Ag  may  be  expressed  as 


=  /[/2>''->o"G,(r„rpdgg(r_)2xaVr, 


(20) 


O' 


^eWs  variables  of  integration  have  been  omitted.  Evaluating  Eq. 


(20) 


^  _7tQ^ein(r^/r^) 

®  2p 


(21) 


from  which  Eq.  (19)  yields  the  common  branch  distribution  function. 


4p 


(22) 


It  has  been  deduced  that  the  total  number  of  common  branches 
arbitrary  plane  scales  approximately  with  the  total  number  of  pores  of 
Hence,  there  is  a  probability  of  interconnectivity  at  pore  size  r  that  is 
size.  Defining  this  probability  as  Pi(rp)  via  Eq.  (22),  ^ 


of  size  rp  in  an 
that  size  in  that  plane, 
logarithmic  in  pore 


Pi(r,) 


4p 


(23) 


It  IS  app^ent  that  approximately  one  percent  of  all  pores  of  all  sizes  are  interconnected 
through  larger  pores. 

The  broad  size  range  associated  with  the  interconnectivity  suggests  that  a  very  wide 
range  of  pore  sizes  control  transport  and  that  a  compHcated  mixture  of  convective  and 
diffusive  transport  persists  through  all  of  pore  space.  While  permeability  is  dominated  by  the 
largest  pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  in  smaller 
pores  m  order  to  accurately  descnbe  the  fine  scale  transport  necessary  to  assess  chemical 
reactions.  In  the  following  sections,  bulk  permeability,  bulk  diffusivity  and  small  scale 
convection  are  addressed  using  the  interconnectivity  derived  above. 
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IV.  PERMEABILITY 


The  fundamental  relationship  governing  convection  in  a  porous  media  is  Darcy's  Law 

which  relates  the  volume  flow  rate  to  the  pore  radius  rp  and  the  pressure  gradient  — 

dx 


where  p  is  the  viscosity  of  the  fluid.  The  bulk  permeability  (k)  is  defined  in 
terms  of  the  volume  flow  rate  across  the  cross  sectional  area  A 


Qv^ 

A{-dpldx) 


(25) 


Convection  across  plane  AA  in  Fig.  4  will  p)ossess  contributions  from  two  primary 
sources  illustrated  in  Figs.  4a  and  4b.  Fig.  4a  illustrates  the  case  where  the  convection  in 
plane  AA  is  due  solely  to  the  pores  that  are  interconnected  in  that  plane.  Fig.  4b  illustrates 
the  case  where  the  convection  in  plane  AA  is  due  to  the  smaller  pores  in  the  pore  tree  that 
^e  mterconnect^  outside  of  plane  AA.  This  connectivity  will  translate  into  a  slower  velocity 
in  the  pore  crossing  plane  A  A  but  could  be  significant  because  99%  of  the  pores  in  plane  A  A 
are  not  interconnected  in  that  plane. 


a)  Pores  Interconnected  in  Plane  AA  b)  Pores  Interconnected  Out  of  Plane  AA 

Fig.  4  Convection  in  the  Pore  Tree 
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Consider  any  pore  of  radius  r,  in  plane  A  A  of  Fig.  4b  to  be  the  trunk  of  a  tree.  Each 
pore  of  size  rp  within  the  tree  possesses  the  probability  Pjir^)  of  being  interconnected  and 

each  interconnected  pore  in  the  tree  will  carry  volume  flow  rate  Since  there  are 

N,drp  (Eq.  10;  N,=r//rp‘*)  pores  in  size  range  drp  within  the  tree,  the  total  volume  flow  rate 
through  trunk  r,  in  plane  A  A  becomes 

<?.('•.)  =  /  O/r,)  P,(<>)  dr^  (26) 


or,  to  first  order. 


a(0  =  -  " 


ln(r 


max 


32  pp 


+  H.O.T 


(27) 


Within  this  approximation,  it  is  seen  that  is  identical  to  the  volume  flowing 

through  the  pores  that  are  interconnected  within  plane  AA.  i.e.. 


which  demonstrates  that  all  volume  flow  through  plane  AA  in  pore  size  r^  is  dominated  by  the 
interconnectivity  of  size  r,  in  plane  AA  and  not  by  the  interconnectivity  of  smaller  pores  in 
subsequent  branches  of  the  pore  tree.  Simply  stated:  case  4a  dominates  case  4b. 

Since  all  volume  flow  through  plane  AA  is  limited  by  the  interconnectivity  of  the 
pores  in  that  plane,  Eq.  (24)  for  the  volume  flow  rate  may  be  rewritten  to  include  all 
interconnected  pores  in  area  A.  Hence, 


Q 


(29) 


where  /(r^)  is  the  "common  branch  distribution  function"  given  by  Eq.  (22).  The  bulk 
permeability  (k)  is  then  expressed  as 


k 


— 


(30) 
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Upon  integration,  Eq.  (30)  becomes 


k  = 


I  16p  } 


(31) 


Equation  (31)  resembles  a  dozen  other  expressions  available  in  the  literature  wherein  it 
is  concurred  that  the  bulk  permeability  is  dominated  by  the  largest  pores  in  the  media  but  the 
unknown  value  of  that  permeability  is  simply  replaced  by  an  unknown  pore  size  to  the  second 
power.  Since  the  pore  size  distribution  function  will  be  least  accurate  at  the  extreme  end  of 
the  size  range,  i.e.  at  r^^j^  ,  no  claim  can  possibly  be  made  that  the  numerical  constants  in 
Eq.  (31)  are  in  any  way  superior  to  those  derived  elsewhere.  One  important  advantage  of  the 
extended  pore  tree  model  is  that  it  characterizes  the  distribution  of  permeability  in  pore  space, 
a  feature  that  will  be  important  in  describing  fine  scale  contaminant  transport  and  in-situ 
remediation.  A  second  advantage  is  the  ability  to  assess  statistical  errors  in  the  measurement 
of  the  permeability  as  a  function  of  the  measurement  scale  size.  This  exercise  is  also  a  good 
test  of  the  extended  pore  tree  model. 


Consider  a  soil  sample  with  the  following  physical  characteristics: 


Conductivity: 

Permeability: 

Porosity: 

Pore  Aspect  Ratio: 

InCW/l-mm) 

From  Eq.  (31),  it  follows  that 
Tmax  =  300  pm 

and  subsequently  r„i„  =  20  A.  The  size  of  the  smallest  pore  is  not  an  important  parameter  for 
this  application  but  may  be  readily  adjusted  through  a  minor  variation  in  the  value  of  P  (e.g., 
^iniii  of  the  order  of  100  A,  p  =  10).  The  bold  assertion  made  in  applying  this  pore 
structure  model  to  soil  is  that  the  1/rp^  pore  size  distribution  is  valid  between  r„,i„  and 

To  investigate  the  role  of  the  measurement  scale  size  on  permeability,  consider  the 
largest  pore  r„„  contained  in  the  spherical  sample  of  radius  "a"  as  given  by  Eq.  (1).  Since 
fmax  in  Eq.  (31)  for  the  permeability  represents  the  largest  pore  in  the  media,  the 
corresponding  value  of  "a"  is  denoted  a^na^  nnd  represents  the  largest  sample  size  for  which 
the  pore  sizes  will  scale  with  the  dimensions  of  the  sample.  From  Eqs.  (1)  and  (31) 


V  =2  cm/hr 
k  =0.6  Darcy 
0  =50% 

K„=5 
p  =12 


max 


24ii:,p 

04/3 


fk 


(32) 


Each  sphere  of  radius  will  contain  one  pore  of  size  r„^.  A  20  x  20  grid  of  these  spheres 
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will  be  characterized  by  the  dimension  40a^  and  contain  400  pores  of  size  r^.  Each  of 
these  pores  possess  probability  P,(rp)  of  being  interconnected.  Following  Eq.  (23),  P,(r  )  is 
approximately  0.0025  for  rp  sufficiently  close  to  r^.  Hence,  only  one  of  the  400  largest 
pores  in  this  20  x  20  grid  will  be  interconnected  and  the  error  in  the  measurement  of  the 
permeability  will  correspond  to  the  statistical  error  of  100%  associated  with  that  of  a  sample 
Carrying  this  argument  to  a  200  x  200  grid  of  dimension  400a_^  ,  there  will 
be  100  mterconnected  pores  corresponding  to  a  statistical  error  of  10%.  Similarly,  a  arid  of 
scale  4000a„,„  will  reduce  the  error  to  1  % . 

Figure  5  illusteates  the  predicted  permeability  measurement  error  associated  with  the 
so  sample  characterized  above  (a„^  =0.3  cm).  Note  that  the  errors  associated  with  the 
measurement  of  permeability  become  negligible  as  the  measurement  scale  size  approaches 
several  meters.  This  has  been  confirmed  by  the  infiltration  data  of  Shouse  et.  al.  (1994).  The 
measured  value  of  hydraulic  conductivity  asymptotes  to  2  cm/hr  at  measurement  scales 
greater  than  4  meters.  At  smaller  measurement  scales,  the  inferred  measurement  error  is 
calculated  under  the  assumption  that  the  asymptote  is  precisely  2  cm/hr.  The  excellent 
agreement  between  the  predicted  and  inferred  error  supports  the  extension  of  the  pore  tree 
model  to  describe  a  porous  permeable  media. 


Figure  5.  Error  in  Permeability  Associated  with  the  Measurement  Scale 
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V.  BULK  GASEOUS  DIFFUSION  IN  A  PARTIALLY  SATURATED  MEDIA 


The  extended  pore  tree  model  is  readily  adapted  to  a  partially  saturated  media  through 
the  assumption  that  all  of  the  water  is  contained  in  pore  sizes  between  r^  and  r^g,  while  only 
gas  is  contained  between  r^^j  and  r^.  Since,  by  Eq.  (6),  porosity  is  distributed  as  In(rp)  in 
pore  space,  the  porosity  associated  with  the  air  filled  pores  (9J  is  approximated  by 


0 


(33) 


where  Tja,  is  treated  as  an  independent  variable  of  the  saturated  pore  structure.  No  gaseous 
diffusion  is  allowed  within  rp  <  r^^,.  With  this  restriction,  the  extended  pore  tree  model  is  used 
to  develop  an  explicit  relationship  between  bulk  gaseous  diffusivity  and  the  permeability 
which  is  validated  through  the  diffusivity  data  of  Washington  et.  aJ.,  (1994).  It  is 
demonstrated  that  the  gas  diffusivity  scales  as  l/r^^j  and  it  is  the  sensitivity  of  r„,  to  the 
saturated  volume  that  controls  the  saturated  diffusivity. 

The  diffusive  mass  flux  in  a  single  pore  is  given  by 

where  Dg  is  the  continuum  gas  diffusion  coefficient  (Dg=0.2  cm^/s).  Continuum  gas  phase 
diffusion  is  valid  only  for  pore  radii  sufficiently  large  that  Dg  is  greater  than  the  Knudsen 
diffusion  coefficient,  DKn=2Vrp/3,  where  V  is  the  mean  thermal  speed  of  a  molecule. 

Knudsen  diffusion  is,  characterized  by  gas  collisions  with  the  jwre  walls  and  is  valid  only  for 

rp  <  3Dg/2V  =  0.1|j,m.  The  limit  of  validity  of  Eq.  (34)  is  denoted  by  ,  the  larger  of 
either  r^^,  or  3Dg/2V.  This  limit  restricts  continuum  gaseous  diffusion  from  both  the  saturated 
pores  and  the  unsaturated  pores  controlled  by  free  molecule  flow. 

Just  as  in  the  case  of  convection,  it  must  be  determined  whether  the  mass  flux  across 
plane  AA  in  Fig.  4  is  determined  by  the  interconnectivity  of  the  smaller  pores  out  of  the 
plane  (Fig.  4b)  or  by  only  those  pores  that  are  interconnected  in  the  plane  (Fig.  4a).  Consider 
any  pore  of  radius  r^  in  plane  AA  of  Fig.  4b  to  be  the  trunk  of  a  tree.  Each  pore  of  size  rp 
within  the  tree  possesses  the  probability  of  being  interconnected  and  each 

interconnected  pore  in  the  tree  will  carry  the  mass  flow  rate  Mp(r^)  .  Since  there  are  N^drp 

(Eq.  10:  N,=rsVrp‘‘)  pores  in  size  range  drp  within  the  tree,  the  total  mass  flow  rate 
through  each  and  every  trunk  of  radius  r^  in  plane  AA  becomes 

M.(r,)  .  /  M/r,)  P,(rp  N,  dr^  (35) 
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Integration  of  Eq.  (35)  yields  the  mass  flux  (case  4b)  for  each  trunk  of  radius  r^ 


f-ap< 


4r 


dx 


sat 


(36) 


where  a  ln(r^^/r,J  term  was  eliminated  via  Eq.  (33). 

If  the  mass  flux  through  plane  AA  is  limited  by  the  pores  that  are  interconnected  in 

that  plane  (case  4a),  the  mass  flux  is  expressed  as  Pj(r^) ,  and  it  is  immediately  seen 

that 


(37) 


i.e. ,  case  4b  dominates  case  4a.  Since  the  mass  diffusion  in  plane  AA  is  determined  by  the 
pore  interconnectivity  of  the  smaller  pores  outside  of  plane  AA,  the  saturation  of  those 
smaller  pores  becomes  an  important  element  in  the  bulk  gaseous  diffusion. 

Since  each  pore  of  radius  r,  in  plane  A  A  carries  mass  flux  the  bulk  diffusion 

coefficient  is  obtained  by  integrating  Eq.  (36)  over  all  pores  in  that  plane 


or,  upon  integration 


D 


bulk 


7  _ 

■  f 


4r 


sat 


D  00  r 

n  -  8  a 

^bulk 


sat 


(38) 


(39) 


and  upon  eliminating  r^^,;  via  Eq.  (31),  the  bulk  diffusivity  is  expressed  in  terms  of  the 
permeability  k. 


D 


bulk 


'  sat 


(40) 


The  bulk  diffusivity  cannot  increase  indefinitely  with  increasing  permeability  as 
inferred  by  Eq.  (40).  In  deriving  this  expression,  the  mass  flux  through  the  interconnected 
branches  of  the  tree  was  not  constrained  from  exceeding  the  diffusive  capabilities  of  the  trunk 


14 


itself.  To  correct  this  potential  problem,  the  limit  of  is  determined  as  the  maximum 
diffusive  flux  in  pore  r,  integrated  over  all  pores  in  plane  AA. 


^xoku 

This  limit  is  illustrated  in  Figure  6  together  with  the  predicted  values  of  0^^^,  for  an 
extended  range  of  values  of  permeability  and  the  saturation  radius,  r,^.  Model  predictions 
correspond  to  the  measured  values  of  6, =0.2  and  6=0.5  from  Washington  et.  al.,  (1994), 
and  the  diffusivity  data  suggest  a  value  of  r,„  in  the  range  of  10  pm  to  100  pm.  An  exact 
comparison  of  the  present  theory  to  the  least  squared  fit  of  the  data  suggests  a  value  of 
30  pm.  While  the  excellent  agreement  with  the  data  of  Washington  et.  al.,  (1994)  does 
substantiate  the  present  theory,  there  is  clearly  a  very  wide  range  of  possible  values  for 
in  partially  saturated  soil  which  will  depend  upon  an  unknown  saturation  radius.  A  two  order 
of  magnitude  decrease  in  the  saturation  radius  will  increase  the  bulk  diffusivity  by  two  orders 
of  magnitude  and  yet  the  corresponding  increase  in  the  air  filled  porosity  is,  by  (33),  only 
33%.  Hence,  field  measurements  of  the  unsaturated  volume  are  not  sufficiently  accurate  to 
correlate  bulk  diffusivities.  If  bulk  diffusivities  are  to  be  correlated  with  field  data,  such 
measurements  should  attempt  to  measure  the  saturation  radius. 


10”  10”  10'*  10^  10^  10^  lO'*  10*  10"* 


PERMEABILITY  (cm*) 

Figure  6.  Bulk  Diffusivity  of  a  Partially  Saturated  Soil 
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VI.  PERMEABLE  SUB-RANGE 


The  present  theory  offers  the  opportunity  to  investigate  some  of  the  features  of  the 
pore  structure  that  are  relevant  to  the  fine  scale  contaminant  transport  and  remediation 
reactions.  One  such  feature  is  the  distribution  of  permeability  in  pore  space.  While  the  bulk 
permeability  is  dominated  by  the  largest  pores,  some  permeability  occurs  in  smaller  pores.  It 
is  the  balance  of  the  sub-scade  convection  with  the  small  scale  diffusion  that  will  control 
contaminant  transport  and  in-situ  remediation. 

The  analysis  and  data  comparison  on  the  errors  introduced  by  the  measurement  scale 
size,  section  IV,  has  shown  that  soil  with  a  bulk  permeability  of  0.6  Darcy  possesses,  on 
average,  one  large  interconnected  pore  in  any  12  cm  by  12  cm  cross  section.  To  allow  the 
penetration  of  the  permeate  throughout  such  a  coarse  grid,  it  will  be  shown  that  there  is  a 
very  extensive  permeable  sub-range  in  the  pore  structure  which  does  not  contribute 
significantly  to  the  bulk  permeability  but  in  which  convection  dominates  diffusion. 

To  evaluate  the  permeable  sub-range,  Eq.  (31)  is  rewritten  to  express  the  permeability 
contained  in  pores  of  radius  rp  (rp  <  r,„ax)  or  smaller 


K  = 


16  p 


\2 


(42) 


where  k„  represents  the  permeability  in  the  sub-range  dimension  to  be  defined  by  L„(rp).  It 
was  deduced  above  that  L„(rp)  must  be  selected  such  that  there  is  at  least  one  interconnected 
pore  of  radius  rp  in  cross-sectional  area  in  order  for  k^^  to  be  accurate  to  order  unity. 
Subsequently,  a  20  x  20  grid  of  spheres  of  radius  a,,  will  contain  one  interconnected  pore  of 
radius  rp 


L  ( r  )  =  4Qa 

sr^  p' 


where  a„  is  given  by  Eq.  (32)  with  k  replaced  with  k^^. 


(43) 


(44) 


Eqs.  (42)  to  (44)  represent  the  scaling  of  the  pore  structure  with  scale  size  for  all  rp  <  r^,^  or 
for  all  a  <  and  all  Lj^tp)  ^  L„(r„,„)  .  The  balance  of  convection  and  diffusion  within 
these  scale  sizes  may  be  assessed. 

The  convective  flux  of  species  "c"  with  mass  density  p<.  across  area  is  denoted  by 


=  lI  Pc(  -dpldx)l  p 


(45) 
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and  the  diffusive  mass  flux  across  the  same  cross  section  is  expressed  as 


where  the  characteristic  gradient  of  p,  is  pj  L„  and  the  effective  diffusion  coefficient  is 

The  balance  between  convection  and  diffusion  occurs  on  the  length  scale  L  where  the 
convective  and  diffusive  mass  fluxes  are  equal.  From  Eqs.  (45)  and  (46) 


L  =  100 


Q^(,-dpldx) 


1/3 


(47) 


The  boundary  between  convective  and  diffusive  transport  is  illustrated  in  Figure  7  for  a  solute 
in  a  saturated  porous  media.  The  effective  diffusion  coefficient  is  assumed  to  be  given  by  the 
porosity  times  the  solute  self  diffusion  coefficient  D,  (D,~  10’^  cmVs).  The  boundary  between 


Figure  7.  Permeable  Sub-Range  for  Convection  in  a  Saturated  Porous  Media 
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sub-scale  convection  and  diffusion  is  illustrated  for  the  solute  in  water  (|j.  =  0.01  poise)  at 
pressure  gradients  of  0.01  atm/m,  0.1  atm/m  (hydraulic)  and  10  atm/m.  Also  illustrated  is  the 
length  scale  associated  with  the  spacing  of  the  large  pores  controlling  the  bulk  permeability 
(Eqs.  (43)  and  (44)  in  the  limit  of  rp  — >  ).  The  results  depict  a  very  extensive  permeable 

sub-range  at  bulk  permeabilities  characteristic  of  soil,  sand  and  gravel.  For  soil  with  a  bulk 
permeability  of  1  Darcy,  the  large  interconnected  peres  are  nominally  15  cm  apart.  The  soil 
matrix  between  the  large  peres  experiences  a  convective  mass  flux  through  its  smaller  pores 
whose  mass  flow  rates  are  much  smaller  than  that  of  the  bulk  permeability  but  greater  than 
that  possible  by  diffusion.  At  length  scales  of  order  one  centimeter,  diffusion  becomes  rate 
limiting.  Hence,  this  px)re  structure  model  provides  a  methodology  for  determining  the  length 
scale  separating  the  mobile  and  immobile  regions  of  the  soil.  The  theory  depicts  that  the  size 
of  the  immobile  region  is  primarily  a  function  of  the  pressure  gradient. 

When  the  above  analysis  is  applied  to  gas  phase  diffusion  in  an  unsaturated  media, 
pi  is  of  the  order  of  10  '*  px)ise  and  is  0.05  cm^/s  (Figure  6).  The  boundary  between 
convection  and  diffusion  occurs  on  length  scales  five  times  as  large  as  those  depicted  in 
Figure  7  for  the  same  pressure  gradients.  While  the  size  of  the  permeable  sub-range  is 
severely  reduced,  the  entire  distribution  of  sp>ecies  "c"  is  still  limited  by  a  diffusive  length 
scale  that  is  primarily  a  function  of  the  gas  pressure  gradient.  Hence,  the  length  scale 
separating  the  mobile  and  immobile  regions  of  the  soil  is  of  the  order  of  five  centimeters  at 
gas  pressure  gradients  equal  to  the  hydraulic  pressure  gradient  (0.1  atm/m)  and  of  the  order 
of  10  cm  at  0.01  atm/m. 

Similarly,  when  the  above  analysis  is  applied  to  gas  phase  diffusion  in  a  partially 
saturated  media,  pi  is  of  order  10-“  poise  and  may  be  as  low  as  lO'^  cmVs  (Figure  6). 
Under  these  conditions,  the  boundary  between  convection  and  diffusion  occurs  at  length  scales 
identical  to  those  depicted  in  Figure  7  and  the  existence  of  a  p)ermeable  sub-range  is 
preserved.  However,  for  a  partially  saturated  media,  Db„u,  should  be  used  for  .  From  Eq. 
(40)  we  write 


(48) 


from  which  it  follows  that  the  length  scale  separating  the  mobile  and  immobile  regions  of  a 
partially  saturated  soil  has  a  slightly  stronger  pressure  gradient  dejjendence 


L 


=  44 


e^a-dp/cb) 


V2 


(49) 


It  is  anticipated  that  the  concept  of  the  permeable  sub-range  will  help  researchers 
develop  relatively  simple,  physics  based  submodels  for  the  ground  water/remediation  codes.  If 
these  submodels  were  tested  independently  from  the  codes,  the  codes  themselves  would  not 
require  parameter  "fitting"  and  could  become  more  directive  than  interactive. 


18 


VII.  SUMMARY 

The  pore  tree  model  has  been 

characterizes  the  subsurface  S  sites  The  random  nature  of  the  pore  structure, 

and  the  in-situ  remediation  of  ^^^tamin^  to  porous  sod 

which  formed  the  basis  of  the  stabstica  •  obtained  via  a  statistical  det^mmation 

and  sand.  The  interconnectivity  of  the  convection  and  bulk  diffusion 

of  the  "branches"  that  are  common  to  seve  diffusion  and  coupled  chemical 

through  the  large  scale  (mobile)  structure.  The  statistical  analysis  reported  above 

reactions  within  the  smaller  sc^^^  across  the  entire  pore 

has  determined  that  the  probability  accompanying  a  decreasing  pore  size, 

size  range,  with  a  slight  increase  in  it  if  also  important  to  establish  the  level 

While  permeability  is  dominated  y  |  intermediate  scales  in  order  to  accurately 

seal!  difhision  and  coupled  chemical  reactions. 

The  permeability  across  a 

interconnected  in  that  plane.  The  T^J^^on^cted  This  establishes  a  very  coarse 

quarter  of  one  percent  of  all  large  ^res  si^e  errors.  The  extend^  pore 

grid  for  the  permeability  which  leads  g^jent  errors  in  the  permeability  of  sod  due 

free  model  has  successfully  ^9^^  has  indirectly  confirmed  the  low 

to  the  measurement  scale  size  (Shouse,  et.ai.,  i:^:^  ; 
probability  of  the  interconnectivity. 

•  it.,  nrross  a  eiven  plane  is  shown  to  be  limited  by  the 
The  bulk  gaseous  diffusivity  across  a  giv  p  be 

toterconnectivity  of  the  ■>"  “'r^foQ4t 

saturated,  resulting  in  a  strong  diffusivity  data  of  Washington  et.  al.,  (1994) 

pore.  A  comparison  of  the  preset  agreement  with  the  data  does 

suggests  a  saturation  radius  of  30  0™- ^  pj^iaily  saturated  soil  is  very  sensitive  to  an 

substantiate  the  P«*”' diSsivifcs  are  to  J  correlated  with  field  data,  such 
— mS™Cldt;:m‘[rr"e  the  saturation  radius. 

The  permeabUity  and  the  bulk 

structure  and  pore  interconnectivity  by  the  out  of  plane  (Fig.  4b) 

Sonnectiv^  (Fig.4a)  and  by  P^r^^nectivity  and  bulk  diftision 

interconnecbvity.  Permeability  is  i  The  apparent  success  of  th^ 

is  limited  by  the  “''’‘5’ that  the  extended  pore  tree 

rSlyTsc^r^- — " 

transport. 

Analysis  of  the  permeability  ofT 

the  distribution  of  the  hich  does  not  contribute  significantly  to  the  bulk 

permeable  sub-range  in  the  P°^®  •  ^tes  diffusion.  It  is  the  balance  of  the  sub-scale 

"l  "Zion  th.  wm  comrol  contaminant  transport  and 
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-bile 

or  size  basis,  moi  reaSite  Z  Z^er  aZf  h' 

counterparts.  Simultaneously  Ae^n!  convective  flow  than  their  larger 

resistance  to  convection  and  a  reduced  arn^^y  surrounding  the  smaller  grains  offer  more 

available  in  the  convective  flow  The  size  remediation  chemicals  will  be 

variations  in  fluid  veS  Snl^d  ciistnbution  of  the  pores  and  grains,  and  the 

these  transport  proc:^TFT:Zs1Z2  Tdevll^ 

convection  and  small  scale  difiusiOT  uoon  which  ^  methodology  to  couple 

accuralely  describe  — cZZls^TriZSt-sr^^ZS^  “  ^  “> 
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